Constraint on the Low Energy Constants of Wilson Chiral Perturbation Theory 
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Wilson chiral perturbation theory (WChPT) is the effective field theory describing the long- 
distance properties of lattice QCD with Wilson or twisted-mass fermions. We consider here WChPT 
for the theory with two light flavors of Wilson fermions or a single light twisted-mass fermion. 
Discretization errors introduce three low energy constants (LECs) into partially quenched WChPT 
at 0{a^), conventionally called Wq, and Wg. The phase structure of the theory at non-zero 
a depends on the sign of the combination 2VKg -I- Wg, while the spectrum of the lattice Hermitian 
Wilson-Dirac operator depends on all three constants. It has been argued, based on the positivity of 
partition functions of fixed topological charge, and on the convergence of graded group integrals that 
arise in the e-regime of ChPT, that there is a constraint on the LECs arising from the underlying 
lattice theory. In particular, for Wl = = 0, the constraint found is VKg < 0. Here we provide 
an alternative line of argument, based on mass inequalities for the underlying partially quenched 
theory. We find that Wg < 0, irrespective of the values of Wg and Wj. Our constraint implies that 
2Wfi > |Wg| if the phase diagram is to be described by the first-order scenario, as recent simulations 
suggest is the case for some choices of action. 
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I. INTRODUCTION 

Effective field theories such as chiral perturbation the- 
ory (ChPT) contain coefficients, usually called low energy 
constants (LECs), which are not determined by symme- 
try. If the matching between the high and low-energy 
theories is non-perturbative, as is the case in the match- 
ing of QCD to ChPT, then the LECs must be determined 
either by experiment or by a non-perturbative method 
such as lattice QCD. One usually has no information 
on the LECs, other than a prediction for their order of 
magnitude based on naive dimensional analysis. It is 
sometimes possible, however, to constrain the signs of 
particular LECs based on the physics of the high-energy 
theory. For example, certain four-derivative terms in the 
chiral Lagrangian are constrained to be positive based on 
causality [l| . This argument has been generalized and ap- 
plied widely in Ref. [2|. Another example concerns the 
chiral Lagrangian describing a lattice simulation at non- 
zero lattice spacing with a mixed-action (different valence 
and sea-quark actions). It is found in Ref. [3] that, using 
generalized QCD mass inequalities one finds a con- 
straint on a combination of the LECs which arise due to 
discretization errors. 

A further method of constraining LECs has recently 
been discovered in the context of calculating the low- 
energy spectrum and eigenvalue properties of the lat- 
tice Hermitian Wilson-Dirac operator [5, ji]. One line 
of argument is based on the positivity of the underly- 
ing two- flavor fermion determinant, which follows from 



the 75— hermiticity of the Wilson-Dirac operator. Specif- 
ically, the partition function at fixed (odd) topology is 
positive in the underlying theory but is only positive 
in the effective theory (here partially quenched Wilson 
ChPT [PQWChPT]) if the LECs satisfy a constraint Q. 
In the standard convention for LECs,""^ this constraint 
is < if ^ ^ 0.2 Another line of argu- 
ment notes that the partially-quenched partition func- 
tion for zero-momentum modes (which determines the 
leading order behavior in the e-regime) converges only if 
< 1^6 + H-^ We also note that a similar con- 
straint (specifically, < independent of Wg and W^) 
was found earlier by one of us, based on the finding that 
the method for calculating the spectral density in infinite 
volume using PQWChPT only worked if the constraint 
held 0. It was not clear, however, whether this was a 
fundamental constraint or simply a shortcoming of the 
method of calculation. 

The constraints found in Refs. [1, imply an in- 
teresting corollary if one assumes that scaling at large 
Nc (number of colors) is a good guide at Nc = 3. In 
particular, since M^g/^s ~ ^j/Wg ^ l/Nc, this as- 
sumption would mean that one can ignore Wq and W^^ 
to first approximation. Then the constraints imply that 
any discretization of Wilson fermions will have an Aoki 
phase for small enough physical quark mass. The other 
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^ Note that our convention for W^, which follows Ref. Jvj, differs 

in sign from that used in Refs. 0,01. 
^ Since these LECs appear at leading order in the appropriate 

power-counting, they are independent of the renormalization 

scale. 

^ It may be possible to obtain further constraints from these or 
similar lines of argument @, 0. 
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possible phase diagram — the first-order scenario [lOj — 
would not occur. This is in apparent contradiction with 
the results of simulations using dynamical twisted-mass 
fermions, which find strong evidence for the first-order 
scenario |llH20ij . Of course, large Nc scaling may not be 
useful for Nc = 3, in which case the connection between 
the constraint and the phase scenario need not hold. 

Given this situation, we think that it is important 
to find an alternative line of argument leading to such 
constraints. This is what we provide in the present 
note. In particular, we find that a generalization of 
the mass-inequality method of Ref. [3| constrains the 
LECs of WChPT.^ Our constraint results from consid- 
ering the twisted-mass generalization of Wilson fermions 
and comparing the quark- connected part of the neutral 
pion propagator to the charged pion propagator. A par- 
tially quenched set-up is required to separate the quark- 
connected and disconnected contractions, and this is why 
it is the LECs of partially quenched WChPT which enter. 
We find W^8 < 0, independent of and W^. 

The remainder of this note is organized as follows. In 
the following section we explain how partial quenching al- 
lows one to separately calculate the quark-connected part 
of the neutral pion correlator. In Sec. IIIII we present the 
calculation of the quark-connected neutral pion "mass" 
at leading order in WChPT. We do so only at maximal 
twist, since this suffices to show the constraint. In Sec. lIVI 
we derive an inequality among quark-connected correla- 
tion functions, from which follows the above-noted con- 
straint. We summarize and offer some concluding com- 
ments in Sec. El We relegate some technical details to 
two appendices, the first concerning the form of the con- 
densate in the partially quenched theory, and the second 
extending the analysis of the main text from maximal to 
arbitrary twist. 



II. USING PARTIAL QUENCHING TO SELECT 
QUARK-CONNECTED CORRELATORS 

Our argument uses twisted-mass fermions [12, [2^ , so 
we begin by recalling the salient features of this approach. 
In an unquenched theory, the quark Lagrangian takes the 
form 

Cq = qs{Dw + mo + ifio'^5T3)qs , (1) 

where qs is an isodoublet of quark fields (corresponding 
to the up and down quarks), and Dw is the Wilson- 
Dirac operator. The subscript "S" indicates that these 
are sea-quarks, appearing in the fermion determinant, as 
opposed to the valence quarks introduced below. ^ We 



* The fact that mass-inequalities can provide useful information in 
twisted-mass theories has also been noted in Refs. [sl, |21|| . 

^ Simulations using two doublets of dynamical twisted-mass 
fermions are also now being done, with the second such fermion 



refer to mo as the normal (bare) mass and fiQ as the 
twisted (bare) mass. The following considerations do not 
depend on whether Dw is improved, or on the form of 
the gauge action, so we do not specify either. We will 
need only the property of "75-hermiticity" : 

l5Dwl5 = . (2) 

When writing the Lagrangian in the form ([T}, we are 
using what is commonly called the "twisted basis", in 
which the mass, and not the Wilson term, is twisted. 

In the continuum limit, a mass term 771-1-1/175X3 can be 
rotated into a purely normal mass ruq = \J vn? + by an 
appropriate axial rotation. Thus the apparent breaking 
of flavor by the /i term is misleading — flavor is preserved 
for all fi. At non-zero lattice spacing, however, flavor 
is explicitly broken from SU{2) to C/(l), leading to a 
splitting of the pion multiplet: 77i^± ^ ttIttQ. As is well 
known (and as will be seen explicitly in the following 
section) the splitting is of 0{a?). 

The particular case of maximal twist corresponds to 
tuning 777o — )■ 777c such that the physical normal mass 
vanishes (or is, at least, sufficiently small compared to 
the twisted mass) . There are a number of different tuning 
criteria that can be used, leading to results for physical 
quantities differing only at 0(af_). For discussion of these 
issues see Refs. [TtI. ITsI. [20I [25l43l| . All that matters here, 
however, is that a consistent criterion exists in which 777^ 
is fixed, such as the one based on the PCAC mass used 
in practice in present simulations [T8l - [20| . 

We will also need to know the quark-level operators 
which couple to the charged and neutral pions in the 
twisted basis. These are given, e.g., in Appendix A of 
Ref. |18i] . The charged pions are created by 

= iqsl5T±qs [t± = -^(ri ± 7T2)] (3) 

(independent of twist angle), while the neutral pion is 
created at maximal twist by 

^" = -qsqs ■ (4) 
Thus the two-point correlators of the charged fields, 

C±(77) = (PT(0)P±(77)) (5) 

(77 labeling lattice sites) , have only quark-connected con- 
tributions. For example, 

C+(77) = 2(tr (75G(m)o,«75G(-m)„,o)) , (6) 



describing the strange and charm quarks [2Cll | . The arguments in 
this note apply equally well to such a set-up, however, because 
the second doublet contains degrees of freedom that are heavy 
on the scale of the light up and down quarks. Thus the form of 
the chiral Lagrangian used in Scc. lIIIl is unchanged (although the 
values of the LECs will be different), and the argument for the 
mass inequalities in Sec. lIVI goes through unchanged. It is impor- 
tant in this regard that the determinant in such Nf = 2 + 1 + 1 
simulations remains real and positive l2i . 
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where the trace is over (imphcit) color and Dirac indices, 
and the quark propagator is 

The neutral pion propagator, however, has both quark- 
connected and disconnected contributions 

= (5°(0)5°(n)) (8) 

+G{-^i)o.,M~^^)n,o)) (10) 

= (tr(GO.)o,o + G(-/i)o,o) 

xtr(G(M)„,„ + G(-/x)„,„)).(ll) 

In practice, the quark-connected part has a much better 
signal to noise ratio than the disconnected part, but im- 
proved techniques have allowed the computation of the 
latter with errors small enough that the mass of the neu- 
tral pion can be extracted [18l |. 

What we are interested in here, however, is the quark- 
connected part of the correlator. Since this is measured 
with small errors (comparable to those for the charged 
correlator) it is worthwhile investigating whether it con- 
tains useful information. In the physical two-flavor the- 
ory one cannot separate the two Wick contractions. It 
is well known, however, that if one considers the par- 
tially quenched (PQ) extension of the theory [s^, then, 
by adding enough valence quarks, one can find correla- 
tion functions which pick out any desired Wick contrac- 
tion. In the present case it suffices to add a single valence 
isodoublet qv and the corresponding ghost quark isodou- 
blet qy . The Lagrangian for each of these quarks is the 
same as that for qs [the twisted-mass Lagrangian ([1])], 
except that qs is replaced by q^ for the ghost quark. ^ 



Within this PQ setup, the correlation function which 
yields the quark-connected part of the neutral pion cor- 
relator involves the mixed valence-sea pion:^ 

G°'™'-(n) = (qsqviO) qvqs{n)) . (12) 

This is because there is no disconnected Wick contraction 
between a qs and qv- Although we will not use it, it is 
perhaps of interest to note that the disconnected neutral 
pion Wick contraction can be obtained as 

GO'd'-(n) = {qsqsiO) qvqv{n)) . (13) 

Partial quenching is often used to consider valence 
masses (or actions) differing from those of the sea quarks. 
In this work, by contrast, the valence and sea quarks 
have identical actions and masses. Thus there is an exact 
SU{2) flavor symmetry mixing valence and sea quarks. 
This is a subgroup of the SU{4:) flavor symmetry that 
emerges in the continuum limit (itself a subgroup of the 
graded flavor group U{4:\2) that holds for perturbative 
calculations in the continuum PQ theory [35|). 

We also remark that G*'''^°"" and C°'<^^^'^ are separately 
unphysical — they cannot be expressed in terms of a sum 
of exponentially falling terms with positive (real) coef- 
ficients. Nevertheless, they can be calculated using the 
appropriate low-energy effective theory — PQWChPT — 
which itself is an unphysical theory, although perfectly 
well defined in Euclidean space. It turns out that G°''^°"" 
does have, at leading order in WChPT, a physical form 
at long distances, which is all that we need for our argu- 
ment. 



III. WILSON CHPT CALCULATION OF 
CONNECTED PION MASSES 

In this section we calculate (7°''=°"'^ and C± usmg 
PQWChPT. We are interested in the differences between 
these two correlators, which turn out to arise at 0{a^). 
Thus we work to leading order (LO) in the "large cut-off 
effects" or "Aoki" regime in which the power counting 
is m ^ /i ^ o^Aqqy), where m and /i are the physical 



This glosses over an important subtlety. In the ghost sector, con- 
vergence of the functional integral requires that the real part of 
the eigenvalues of the discretized fermion operator are positive. 
This is not the case for Dpi^ -I- mo -|-i/^T375 given that one always 
works with mo < 0. This issue has been resolved, in the context 
of the quenched theory, in Ref . |33|| , and a simple generalization 
works here. The solution is to do an axial rotation in the 
direction by angle it /A, such that the fermion operator becomes 
D — i-f5T3{W + mo) + fJ., where D is the naive discretization of 
the Dirac operator and W the Wilson term. This new operator 
consists of an antihermitian part, with purely imaginary eigen- 
values, and a real offset fj,, which we choose to be positive. (For 
negative fi, an axial rotation in the other direction resolves the 
problem.) For maximal twist, this is exactly the axial rotation 
that brings one to the physical basis, but for other twist angles 
it gives a different basis. In this new basis one can add in va- 
lence and ghost quarks. One then goes through the standard 
steps to obtain the chiral Lagrangian including discretization 



errors llOll. following a simple generalization of the analysis of 
Ref. [331 . Compared to the usual chiral Lagrangian, one has ad- 
ditional factors of ztirs in terms containing spurious coming from 
discretization errors. Thus the Lagrangian looks non-standard. 
In the quark sector (sea and valence) one can, however, undo 
the axial rotation (now at the level of the chiral fields) ending 
up with the standard form of the chiral Lagrangian for WChPT 
[presented below in Eq. I lf4l |]. This does not work in the ghost 
sector, since one is not allowed to do normal axial rotations there. 
This restriction does not, however, effect the present calculation, 
since we only consider correlation functions in the quark sector. 
In fact, the correct procedure in the ghost sector has been worked 
out in Ref. [34| , generalizing the methodology of Ref. [sSj . 
We could just as well add two isodoublets of valence quarks 
(and corresponding ghosts) and use {qvi'lV2{0)qv2qvi{iT-)) ■ The 
choice made in the text is, however, the minimal one. 
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normal and twisted masses [defined in Eq. (fT6)) below]. 
There is no need to work to higher order, since for the 
purposes of constraining the LECs we can imagine that 
m, jj, and are arbitrarily small. 

At leading order, and after shifting the quark mass 
to remove an (a) term, the partially quenched chiral 
Lagrangian is [7l.llO| 

£^ = ^Str(9^S]a^I]t)-:^Str(xtS + Stx) 

-a^W^ [Str (E + St)]' - a'^W^ [Str (E - S^)]' 
-a^W^^Str (E^ + [i]t]2) (14) 

Here E G SU{4\2), "Str" stands for supertrace, x — 
2BqM with M the mass matrix, and a — 2Woa. Bq and / 
are continuum LECs (with the convention 93 MeV), 
while Wo, Wq, W-^ and Wg are LECs associated with dis- 
cretization errors. Since our set-up requires just two va- 
lence quarks and two ghosts, the graded chiral symmetry 
is® SUi^L X SU{A\2)r. 

The mass matrix in is related to the bare masses 
in the underlying quark Lagrangian ([1]). For the un- 
quenched theory, we have 



(15) 



Zg^{mo - mc)/a, fi Zp^fio/a . (16) 



Here LOm is the "input" twist angle, and iriq = ^/ri^?~{~Ji^ 
is the physical mass in the continuum limit. Maximal 
twist corresponds to m = 0. Note that in this case one 
does not need to know the renormalization factors Zs 
and Zp in order to determine the twist angle. For the 
partially quenched theory, the mass matrix, which has 
dimension 6 x 6, is block diagonal, with each block con- 
taining M. As noted in the previous section, this mass 
matrix leaves an unbroken SU{2) symmetry between sea 
and valence quarks. 

We must first determine the orientation of the vacuum. 
So = (0|E|0), taking into account the 0{a'^) terms. In 
the unquenched sector this has been done in Refs. [Sf 
m. Writing 



v-\unqu 
^0 



(17) 



one needs in general to solve a quartic [given in Eq. (|B1I) ] 
to determine ujq, and luq — Um is generically of 0(1). For 



The actual symmetry differs from this due to the constraints 
from convergence of ghost integrals. For perturbative calcula- 
tions, such as those we perform here, one can, however, work as 
if the symmetry is as claimed. This was shown for the continuum 
PQ theory in Ref. [35| . and presumably carries over to Wilson 
PQChPT. In fact, all we need in the present calculation are fluc- 
tuations in the quark sector, and here the appropriate symmetry 
is certainly SU{A). For non-perturbative calculations, however, 
such as those done in Refs. [aQ], one must account for the need 
to have convergent integrals in the ghost sector, which leads to 
a different global group. 



the special case of maximal twist, however, the solution 
is simply luq = ujm = ±7r/2, i.e. the input and output 
twist angles are the same. 

For the partially quenched theory, we argue in Ap- 
pendix that 



So = 

























e<p9 



(18) 



in a 2 X 2 block notation with the blocks ordered as sea, 
valence and ghost quarks. In words, this result says that 
the SU{2) valence-sea symmetry is unbroken (one impli- 
cation of which is that the vacuum twist in the valence 
sector is the same as that in the sea sector) and that 
there are no quark-ghost condensates. Wc do not need 
to discuss the (subtle issue) of the ghost condensate e*^" 
(which is a 2 X 2 matrix), since we will not need prop- 
agators involving ghosts. This issue has been discussed, 
albeit in a different power-counting, in Ref. [s^ j. 

Pion masses can now be obtained by considering small 
oscillations around the condensate. To do this, we use 



e*"«^3/2 
e^"°^^/2 
e'*''/2 



(19) 



with Sph = exp(i\/27r//) containing the pion fields. We 
will only need the quark part of the pion field, which we 
decompose as follows 



VqirVq 



T^ss T^sv 

TTVS T^VV 





(20) 



with Vq the projector onto the quark subspace 



Va 




(21) 



The block pion fields in ([20| contain isosinglet compo- 
nents (i.e. ?7-like fields) as well as the usual isovector pi- 
ous, but the isosinglet parts play no role in the following 
calculation. As we show below (following Ref. [s^), the 
symmetric positioning of the condensate in ()19|) leads to 
the usual identification of the individual pion fields in tt. 
In particular, for the isovector fields, the decomposition 
for each block is the usual one, 



7r7V2 7r+ 



(22) 



To show this, we next need to map the operators P^, 
5° and qsqv of Eqs. ([3]), (g]) and ([HI) into the effective 
theory. This is a standard exercise, requiring the intro- 
duction of scalar and pseudoscalar sources into the mass 
matrix M . At LO in our power counting, the results are 
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the same as in the continuum. In particular for quark 
bihnears we have 



Str (PgTpt - S]) (23) 



and 



qTq 



Str(-PgT[I] + I]t] 



(24) 



where T is an arbitrary flavor matrix acting on the 4x4 
quark subspace. Using this result, and the expansion 
(HH), we find 



Str(7'^Sr±[4j^~Sph]) (25) 

-2/i3o7r|5 + 0(^3) (26) 

^Str(7'^^{c„[I]p, + Et,J 

-*sor3[4h-Sph]}) (27) 

-so2/So7rg5 + 0(^2) (28) 

^Str(7'^^{co[S + I]t] 

-zsoT3[4h-Sph]}) (29) 

-so2/So7rO,s + 0(^') ■ (30) 



where cq = coswq and sq — sinwo- V^^ is the projec- 
tor onto the sea-sea block, and V^^ picks out the off- 
diagonal valence-sea block: 



-qsqv 



^SS 














1 




: 





• 




• 





• 



















(31) 



Thus, at maximal twist (sq — 1), and 5° indeed 
couple with equal strength to the charged and neutral 
pions, respectively, as required by the underlying theory. 
We also see that iTyg is the appropriate field to use to 
determine the connected part of the neutral correlator. 

We can now calculate the correlators C^(x), C'^{x) 
and C"'^™""(x) of Eqs. ©, dU and (HI]), respectively.^ 
At this stage we specialize to maximal twist. This not 
only simplifies the resulting expressions but also turns 
out, as sketched in Appendix [B] to give the same con- 
straint on the LECs as one finds when working at arbi- 
trary twist. Expressed in terms of the rotated fields, the 



^ Note that we are now in a continuum theory, so the lattice label 
n is replaced by Euclidean position x (with the correspondence 
X ~ an). 



chiral Lagrangian becomes 
P 



C 



Str (s^Spha^st^ 



^2BoTOqStr fHph + Spi^ 



Str (EoSph + Sjs|,h 
Str (^SoSph 



ytyt 



a^WgStr ( EoSphSoS 



(32) 



Keeping only terms quadratic in the pion fields, we find 
that the term gives 

- w'qP {-k^ss + '^vv + ghost-terms)^ , (33) 
the Wj term vanishes, and the term becomes 

1. n 



'vv\ 



■ ghost— terms J 
(34) 



Here we are using rescaled, dimensionless LECs 



P 



(35) 



We see that, while the Wg term contributes to the masses 
of all the neutral pions, the Wq term contributes only to 
the neutral particles in the diagonal blocks (and thus not 
to the TTgy mass). This is because of the double-strace 
form of Wq , which means that it gives "hairpin vertices" 
in the usual PQChPT parlance. 

Putting this all together, we find that, at leading order, 
each correlator of interest is proportional to the propa- 
gator of the corresponding pion. In momentum space we 
have (still at maximal twist) 



ApBl 



with i = ±, 0, and "0, conn", where^° 



(36) 



{mlY = {m%sY 

(^0,co„n)2 ^ (^0^)2 



2BoA* (37) 
2Bo^l - {2w'^ + w'^)p (38) 
2SoAt - w'J^ . (39) 



The results for m^± and m^o agree with those of 
Refs. [36l - [38| . while that for the connected neutral pion 
is new. It is the latter result which provides the key 
constraint, as we now explain. 



There is one subtlety in the calculation. As can be seen from 
Eq. I I33I I. there are off-diagonal terms proportional to «)g con- 
necting vr^g to T^YY and ghost terms. These do not contribute, 
however, due to a cancellation between valence and ghost contri- 
butions, as must be the case because, for a purely sea-quark pion, 
we can do the calculation solely in the unquenched WChPT, 
leading to the result stated. 
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IV. MASS INEQUALITY AND THE 
CONSTRAINT ON LEGS 



We begin by rewriting the charged correlators using 
75G'(— ^)75 = G'(/i)t (which follows from 75-hermiticity) 

C+{n) - 2(tr(G(M)o,„G(M)i^o)> (40) 
C-{n) = 2(tr(G(M)S^„G(M)„.o)> . (41) 

These two correlators are equal by charge conjugation 
symmetry, i.e. after averaging over each gauge field and 
its complex conjugate. Note that both correlators are a 
sum over positive definite terms, which leads us to expect 
that they are larger than all other correlators (assuming 
appropriate overall normalization factors). This is the 
basis for the mass-inequality method. 

In thepresent case, we can adapt the argument given 
in Ref. [3i]. We start by noting that, on each gauge con- 
figuration. 



or, equivalently. 



< 



[G(m) 



Oa.nb 



2b,Qa 



(42) 
.(43) 



where a and b are color-Dirac indices. Multiplying out, 
summing over the color-Dirac indices, averaging over con- 
figurations [allowed since the quark determinant is real 
and positive], and using Eqs. ^ and ([TU]). we arrive at 
the key inequality 



C+{n) + C~{n) 



C+{n) > 



\n). 



(44) 



which holds for all n. 

Now, for long distances, we can use the forms predicted 
by PQWChPT, which we know from the previous section 
to be (after Fourier transforming) 

C+{n) cx (m+s)i/2(an)~3/2g-m+san (45) 

C'''™""(n) oc {wPsyf''^{an)-^''^e-<^''''', (46) 

with a common coefficient of proportionality. We stress 
that, although (7°'™"" is unphysical, PQWChPT predicts 
that it has a single-particle exponential fall-off at long 
distances. The only way that pS)) and can be con- 
sistent with the inequality (|44p for n large enough that 
the exponential damping dominates is if 



(47) 



The correlators are real and positive, while (70,00111 jg priori 
only known to be real but of indeterminate sign. The PQWChPT 
result II36I I shows, however, that at long distances cOjCont jg ^^^^ 
positive. Thus we chose to consider the sum G(/x) -|- G{fj,)^ in 
Eq. I I42II . so that C'0.c°nn ^ould appear with a positive sign on the 
right-hand-side of the inequality II44I I. We note for completeness, 
however, that we could also have considered the difference in 
G(ai) - G(^)t in Eq. gg, from which one would deduce that, in 
general, C+{n) > |CO'™""(n)|. 



(48) 



Combining this with the results for the masses from 
PQWChPT, Eqs. dSlI) and dMl), we find that 



(49) 



< W!,<0. 



As shown in Appendix [Bl one finds no other constraints 
on the LECs if one repeats the argument at non-maximal 
twist. 

The inequality (|48| can be directly tested in lattice 
simulations, and present results (see, e.g.. Fig. 6 of 
Ref. J_8J) clearly satisfy the inequality. 

We close this section by noting a relationship between 
the mass inequality (|48| and the analysis of the conden- 
sate given in Appendix [X] One of the conclusions from 
the appendix is that the sea-valence SU(2) symmetry 
cannot be spontaneously broken for non-zero /zq. This 
is consistent with the mass inequality because, if there 
were a mixed sea-valence condensate, then one would ex- 
pect that fluctuations in the sea-valence direction would 
diverge, and thus that (m^y)^ would pass through zero 
and become negative. The mass inequality says that this 
cannot happen while m'gg is positive, as it is expected to 
be for any non-zero fio- 



V. CONCLUSIONS 

We have shown that the sign of one of the LECs in- 
duced in Wilson ChPT by discretization errors can be de- 
termined by combining partially quenched WChPT with 
mass inequalities. The core of the argument is technically 
very simple, requiring only a tree-level computation and 
a simple inequality. The only connection between our ar- 
gument and those given in Refs. is that both require 
the positivity of the determinant. 

The constraint wc find is that < 0, independent of 
the values of and Wj. We find no constraints on the 
latter two LECs. These results are the same as found 
in Ref. [9| , based on the failure of a method to calculate 
the spectral density of the hermitian Wilson-Dirac oper- 
ator. Our constraint is also consistent with that given in 
Ref. 0] based on the positivity of the partition function 
in odd topological sectors {W^ < H = = 0). 
It differs from that found using the convergence of the 
zero- mode partition function, namely Wg < Wg + 0] . 
Whether our result is stronger or weaker than this con- 
straint depends on the signs of Wq and W^. 

We stress that all arguments leading to constraints rely 
on the applicability of partially quenched ChPT. In our 
case, we work in the "p-regime" — i.e. large volumes, with 
only small perturbations around the ground state — while 
Refs. 0, Q work in the e-regime in which the zero- modes 
must be integrated over the entire group manifold. 

Our calculation also provides a simple way of deter- 
mining using the result (valid at maximal twist, and 
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generalized to arbitrary twist in Appendix [B|) 

= j^+0{a\a^ml){5l) 

It appears from recent simulations with twisted-mass 
fermions (see, e.g., Refs. [H, Hsj) that this should give 
a fairly accurate determination. The only concern is 
whether the LO contribution will dominate. It would 
thus be interesting to extend the one-loop calculation 
of Refs. [13, im to the partially quenched theory. It 
would also be interesting to compare results obtained us- 
ing (ICTI) with those from other recently proposed methods 
for determining W^, which are based on using a mixed 
action |42| or on partially quenched pion scattering am- 
plitudes [43| . 

We return now to the implications for the phase struc- 
ture of unquenched twisted-mass fermions. As noted in 
the introduction, this depends on the sign of the combi- 
nation of LECs, 2Wq + Wg. If this combination is nega- 
tive, then one is in the Aoki-phase scenario, which means 
that < m° as long as /xq 7^ [as can been seen from 
Eqs. and The resuhs of Refs. how- 

ever, favor the first-order scenario, with > for 
^0 7^ 0. This means that 2Wq + Wg > 0, which, com- 
bined with the inequality < 0, implies in turn that 
2PFg > \W^\. There is nothing theoretically inconsistent 
with this possibility, but it is somewhat surprising given 
that W^/W^ cx 1/iVc for large N^. 

A related implication of the presence of the first-order 
scenario is that quark-disconnected contributions play an 
important role. It is these contributions which, despite 
being suppressed by 1/Nc, lower the neutral pion mass 
below that of the charged pion. This violation of large Nc 
counting (Zweig's rule) is superficially analogous to the 
situation with the 77' in QCD. Here, however, the effect 
has the opposite sign,^^ and is of O(a^) rather than a 
physical effect. 
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Appendix A: Form of the partially quenched 
condensate 

In this appendix we present the arguments for the form 
of the condensate given in Eq. (|18p . The discussion is 
carried out in the underlying theory. 

We first note that we know from a general argument 
given in Ref. ^ that quark-ghost condensates vanish. 
This leads to the zero entries in the rightmost column 
and bottom row (aside from the bottom-right block). 

Secondly, we show that the valence-sea SU{2) symme- 
try is unbroken, leading to the other zeros in (jlSp . as well 
as the result that the condensate in the valence-valence 
block is the same as that in the sea-sea block. The argu- 
ment is a generalization of the Vafa-Witten theorem on 
the absence of fiavor breaking (45j . A similar argument 
was made in Ref. 0] concerning the absence of flavor- 
breaking in the valence sector alone, but this was depen- 
dent on the fact that the valence sector contained quarks 
with an exact chiral symmetry, so that the Dirac opera- 
tor has a continuum-likc spectrum. In the present case 
we have valence and sea Wilson fermions, with no chiral 
symmetry, so the argumentation is different. In fact, it 
is surprising that one can make such an argument at all, 
since we know that the SU (2) symmetry in the sea-sector 
can be spontaneously broken — this is, after all, what hap- 
pens in the Aoki phase. The key difference here is that we 
are working at non- vanishing twisted mass, which avoids 
the appearance of small eigenvalues of Wilson-Dirac op- 
erator. 

We show first that the sea- valence condensate 

{qsuibqvu) (Ai) 

vanishes. The notation here is that, in each 2x2 block, 
we label the two states by u and d. Thus qvu is the 
valence u quark. We work in a volume V ^ at non-zero 
lattice spacing a, and with /io non-zero. We turn on a 
source term 

-Csourcc = A qVul^qSu , (A2) 

chosen to "push" the condensate in a direction such that 
(|A1[) is non-zero. We then take V — ?> 00, followed by 
A — > 0, and find that (jAip vanishes. This implies the 
absence of spontaneous symmetry breaking. 

Explicitly, a simple calculation yields (up to correc- 
tions proportional to A^) 
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This point has been stressed recently in Ref. [44 



V 



{qsul5qvu{n)) = ^(Tr(757; : 75^; ] ^— 

V \ \ Dw + mo + «Aio75 Dw + mo + ^/io75 



1 



1 



V \ \Q + i^iQ Q + ifio 



A / dX p{X) 



1 



I 



(A3) 

(A4) 
(A5) 



Here the traces are over space, Dirac and color indices, 
Q = j5{Dw + rno) is the hermitian Wilson-Dirac opera- 
tor, which has (real) eigenvalues denoted by A, and p{X) 
is the density of eigenvalues per unit volume after aver- 
aging over gauge fields. Note that we expect p(0) to be 
non- vanis hing in general (which gives rise to the Aoki- 
phase l46{ ) but the presence of /zo 7^ shields us from 
the potential singularity at A = 0. Indeed, the coefhcient 
multiplying A is finite for any non-zero a, since the range 
of the integration over A is finite. Thus the sea-valence 
condensate vanishes when A — )• 0. 

Note that to make this argument we need the eigen- 
value density to be well-defined, and for this we need the 
integration over gauge fields to have a positive weight, 
which is the case for twisted-mass fermions. 

Similar arguments show that all the condensates 
{qsjlbWk) vanish, with j and k running independently 
over u and d. Also, by using different twisted masses 
for valence and sea quarks one can show that conden- 
sates {qsjlzlsk — qvjlBQvk) vanish. For the correspond- 
ing scalar condensates, e.g. {qsj^Vk), one ends up with 
expressions such as 



1 



1 



— ■ ( Tr ( 75 



-15- 



(A6) 



Although an eigenvalue decomposition cannot be used 
here, there is no reason to expect that the coefficient of 
A diverges for non-zero a, given the presence of ^0 7^ 0- 
Assuming so, we find that all sea-valence condensates 
vanish. 



Appendix B: Connected pion masses at arbitrary 
twist angle 

In this appendix we give the values of the masses 
(m^)^, (to°)^ and (771°''^°"")^ at arbitrary twist angle. 
For a lattice theory in the large cut-off effects regime, if 
the mass twist angle uJm is not an integer multiple of 7r/2, 
then it will differ from the twist angle in the condensate. 
Wo- As a result, when the leading order chiral Lagrangian 
is expressed in the physical basis (in terms of Sph), the 



form of both the mass term and the O(a^) terms is al- 
tered by the twist. Expanding the leading order chiral 
Lagrangian to ©(tt^), one fixes the relation between ujm 
and Wo by demanding that the linear term vanish. This 
relation is the same as in the unquenched case [36l - [38| : 



2BoAisin(wo - ojm) = -,f{2w'Q + w'g)soCo . 



(Bl) 



This can be used to rewrite the quadratic terms in 
so that they depend only wo and not uJm- One may then 
read off the masses 



(^0,con„)2 



2BoM 

So 

So 
2Bop 

So 



f<4 . 



(B2) 
(B3) 
(B4) 



The results for the unquenched charged and neutral pions 
agree with those of Refs. [37l.[39l|. 

By generalizing the arguments of Section IIVI one can 
show that the connected neutral mass can be no smaller 
than the charged mass. We sketch the generalization 
briefly. The form of the charged correlator, Eq. (|6]), is 
independent of twist. The neutral correlator does, how- 
ever, depend on twist; the operator used to create the 
neutral sea- valence pion becomes 



qsiJbTse ^ qv ■ 



(B5) 



Here w is the twist angle determined in the simulation, 
either from the input masses, or using one of the other 
possible definitions. It will not matter which definition is 
used, since the inequality holds independent of w. Thus 
the connected neutral correlator becomes 



CS'™""(n) = -(tr(i75e^^-^"G(/i)o,„i75e''^«"G(/i)„,o 

+ Z75e-^^="G(/.)S,„z75e-^^="G(A.)i,o){j36) 

Now, using 

[-G(A.)z75e*^^"° +775e-^^^"°G(A.)t]^^„, ' > 0, (B7) 
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and following similar steps as in the main text, one finds 
that 

C°'™""(n) < C+(n). (B8) 

We stress that this inequality holds separately at each 
value of the input bare masses mo and no, and further- 
more, for fixed mo and fio, it holds for any choice of 
w. When we evaluate the correlators in WChPT it is 
most natural to choose uj = ujq, for then the connected 
neutral correlator couples to the sea- valence neutral pion 
with the same strength as the charged correlator does to 

Using other values of oj leads to a weaker inequality. 



the charged pion.-'^^ This means that the WChPT result 
Eq. ([55)1 still holds, except that the masses which appear 
are now those of Eqs. (|B2IIB4p above. 
Putting this all together, it follows that 

Thus, on the one hand, the result < can be demon- 
strated using any non-zero twist angle, but, on the other, 
working at arbitrary twist docs not provide an additional 
constraint on the LECs. 
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